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摘 要 :基于 断裂 力学 理论 ,应 用 有 限 元 法 分 析 带 有 初始 边缘 直 裂 纹 的 二 维 橡胶 板 在 单 轴 拉 伸 下 的 
撕 列 特性。 采用 随机 顺序 吸附 (random sequential adsorption ,RSA) 算法 在 ABAQUS 软件 中 生成 短 纤 
维 随 机 分 布 的 模型 ,研究 网 格 尺 十 和 积分 路 径 对 积分 的 影响 ;计算 复合 材料 在 不 同 拉 伸 载荷 下 的 
积分 值 ;分 析 了 J 积分 (扩展 方向 平行 于 初始 裂纹 方向 ) 随 短 纤维 含量 、 短 纤维 长 度 、 短 纤维 直径 
和 短 纤维 取向 角 的 变化 关系 。 研 究 结 果 表 明 : 网 格 尺寸 和 积分 路 径 对 积分 有 一 定 的 影响 ;J 积分 
随 拉 伸 载荷 的 增 大 而 增 大 ;在 同一 初始 边缘 直列 纹 (ao =6.25 mm) 下 , 短 纤维 填充 橡胶 复合 材料 抗 
撕 列 性 能 随 短 纤维 含量 、 长 度 和 直径 的 增 大 而 提高 , 随 短 纤维 取向 角 的 增 大 而 降低 ; 短 纤维 在 小 角 
度 范围 内 随机 取向 的 抗 撕 裂 性 优 于 大 角度 范围 内 随机 取向 。 
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Abstract: Based on fracture mechanics theory, FEM was used to analyze the tearing characteristics of two- 
dimensional rubber sheet with initial straight edge crack under uniaxial tension. The random sequence ad- 
sorption ( RSA) algorithm was used to generate the random distribution model in ABAQUS. The influence 
of mesh size and integral path on J integral was studied ,and the J integral in different spreading directions 
was calculated. The relation of J integral (the propagation direction is parallel to the initial crack direc- 
tion) with short fiber content ,length , diameter and fiber orientation angle was analyzed. The results show 


that the mesh size and the integral path have an influence on the J integral. The J integral increases as the 
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tensile load increases. Under the same initial straight edge crack (ao =6.25 mm) ,the tear resistance of 


short fiber reinforced rubber composite increased with the increase of short fiber content, length and diame- 


ter, but decreased with the increase of short fiber orientation angle. The tearing resistance of short fibers 


with random orientation in small angle range is better than that in large angle range. 


Key words :short fiber;rubber;, integral; composite materials; crack 


短 纤 维 -橡胶 复合 材料 (short fiber rubber com- 


没有 。 而 短 纤维 的 细 观 结构 对 橡胶 复合 材料 裂纹 扩 


posite ,SFRC) 是 一 种 性 能 优异 、 应 用 广泛 的 高 分 子 
材料 。 短 纤维 特殊 的 形状 系数 使 其 复合 材料 具备 良 
好 的 抗 撕 裂 性 。 然 而 , 短 纤维 复合 材料 的 使 用 仍然 
受到 限制 ,限制 应 用 的 主要 原因 之 一 是 难以 可 靠 地 
预测 其 机 械 性 能 ,尤其 是 手 裂 竺 性。 文献 [1] 和 采用 
虚拟 裂纹 闭合 技术 ( VCCT) 结合 有 限 单 元 法 ( FEM) 
计算 了 短 纤维 橡胶 复合 材料 单 胞 代表 性 体积 单元 
(RVE) 模 型 纤维 端 部 不 同方 向 裂纹 的 应 变 能 释放 
率 。 文 献 [2] 建 立 了 短 纤维 增强 橡胶 复合 材料 在 横 
向 静 力 和 疲劳 拉 伸 载荷 作用 下 的 界面 损伤 计算 模 
型 ,结合 内 聚 力 模型 和 应 力 疲 劳损 伤 模型 ,模拟 损伤 
的 产生 和 扩展 。 文 献 13] 实 现 了 随机 短 云 杉 纤 维 增 
强 聚 丙烯 复合 材料 断裂 行为 的 数值 模拟 ,并 与 数字 
图 像 相关 法 (DIC ) 实验 结果 进行 对 比 。 文 献 [4] 采 
用 有 限 元 分 析 法 研究 SiC 纤维 长 径 比 和 体积 分 数 对 
金属 基 复 合 材料 力学 性 能 的 作用 。 文 献 L5] 利 用 有 
限 元 数值 模拟 方法 计算 了 短 纤维 复合 材料 的 断裂 万 
性 。 文 献 16] 分 析 了 纤维 含量 对 短 玻璃 纤维 复合 材 
料 断 裂 行为 的 影响 。 文 献 [7] 制 备 了 不 同 直径 、 相 
同 长 径 比 的 SFRC 并 通过 实验 方法 比较 了 其 断裂 强 
度 .断裂 形变 等 力学 性 能 。 文 献 [8-9] 应 用 ABAQUS 
软件 对 橡胶 类 双 材 料 界面 裂纹 扩展 进行 非 线性 有 限 
元 分 析 , 发 现 界面 裂纹 扩展 方向 与 两 种 材料 常数 的 
比率 有 关 。 文 献 [10] 发 现 ZrB,-SiC-C, 复 合 材 料 在 
断裂 过 程 中 表现 出 纤维 拔 出 、 纤 维 桥 联 和 和 裂纹 偏转 
的 增 蔬 机 制 。 文 献 [11] 通 过 电镜 发 现 唆 酰 亚 胺 短 
纤维 在 树脂 中 的 无 序 分 布 以 及 其 对 上 下 层 间 的 桥 联 
作用 是 复合 材料 增 韧 的 主要 原因 。 文 献 [12] 通 过 
试验 结合 SEM 手段 验证 了 纤维 的 脱 黏 和 拔 出 对 无 
纺 布 增强 碳纤维 层 合板 分 层 裂纹 起 到 了 阻碍 作用 ， 
从 而 提升 层 间 断裂 箱 性 。 

由 上 述 研 究 现状 可 以 看 出 ,关于 短 纤 维 复合 材 
料 的 断裂 性 能 已 有 较 多 的 研究 ,而 对 短 纤维 -橡胶 复 
合 材 料 破 坏 的 研究 相对 较 少 , 且 多 采用 传统 的 实验 
方法 或 将 其 视 为 各 向 同性 的 均 质 材料 进行 断裂 仿 
真 , 考 虑 复合 材料 非 均 质 特性 的 断裂 数值 研究 几乎 


展 行为 具有 重要 作用 。 因 此 ,本 研究 基于 细 观 力学 
数值 方法 ,以 ABAQUS 软件 为 平台 结合 PYTHON 语 
言 建立 二 维 随机 短 纤 维 橡胶 复合 材料 的 非 均 质 模 
型 ,计算 橡胶 复合 材料 在 不 同 拉 伸 载 荷 下 的 /积分 
值 ,讨论 网 格 尺 寸 和 积分 路 径 对 / 积分 计算 的 影响 ， 
研究 复合 材料 抗 撕 裂 性 随 短 纤维 含量 (V) 长度 
(7 ) 、 直 径 (q) 变化 关系 ,计算 短 纤维 取向 角 在 服从 
数学 均值 为 ,方差 为 o 的 高 斯 分 布下 的 J 积分 值 
用 以 讨论 取向 角 (a) 对 橡胶 复合 材料 抗 撕 裂 性 的 影 
响 , 分 析 了 短 纤 维 在 某 一 角度 范围 内 随机 取向 的 J 
积分 值 随 角度 范围 的 变化 情况 。 


1 理论 依据 


基于 能 量 概念 的 断裂 力学 理论 已 被 证 明 是 分 析 
弹性 体 断裂 的 有 效 方法 。RIVLIN 和 THOMAS 扩展 
了 Griffith 理论 ,成 功 将 应 变 能 释放 率 6G 应 用 于 橡胶 
材料 的 撕 裂 问题 。 后 来 ,RICE 引入 了 J 积分 判 据 作 
为 断裂 准则 ,以 克服 计算 非 线 性 弹性 或 弹 塑性 体 中 
近 端 应 力 场 的 困难 。 从 能 量 的 角度 来 解释 ,J 积分 
为 裂纹 扩展 单位 面积 时 ,系统 所 释放 的 能 量 。 考 虞 
一 个 含有 贯穿 裂纹 的 平面 问题 ,裂纹 表面 无 外 力作 
用 ,围绕 裂纹 尖端 取 一 逆 时 针 路 径 工 , 沿 该 路 径 丁 作 
下 式 积分 便 可 求 得 / 积分 值 。 

一 人 edr - T ds) (1) 


其 中 :w 为 应 变 能 密度 ;三 为 积分 路 径 ;ds 为 弧 长 ;T 
为 张力 矢量 ;u 为 位 移 矢 量 。 

在 ABAQUS 中 采用 域 积分 方法 来 评估 围 线 积 
分 。 已 有 学 者 证 明 J 积分 判 据 可 以 用 于 橡胶 材料 的 
失效 分 析 ” 。 考 虑 如 图 1 所 示 的 逆 时 针 路 径 工 ( 回 
路 始 于 裂纹 下 表面 终于 上 表面 ) ,对 于 橡胶 材料 的 二 
维 问题 ,J 积分 可 表示 为 


= fawn,ds 一 [ (on,+ Tun ) Ou + 
Tr Tr G7 0% 


0 
(Tn 十 on,) 二 ] ds (2) 
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由 于 裂纹 前 缘 的 应 力 场 存 在 奇异 点 ,常规 有 限 
元 网 格 在 求解 裂纹 尖端 处 了 积分 时 难以 实现 。 为 克 
服 常规 单元 在 求解 裂 尖 上 的 弊端 , BARSOUM 和 
HENSHELL 提出 在 二 维 分 析 中 可 以 通过 退化 单元 得 
到 网 格 奇 点 。 如 图 2 中 的 B1 ~ B5 所 示 ,将 边 中 节点 
移 至 靠近 裂 尖 的 1/4 边 长 处 , 即 可 构造 出 +r 的 应 
力 奇 异 点 。 在 本 研究 的 数值 分 析 中 ,裂纹 前 缘 采 用 
1/4 节点 退化 奇异 单元 进行 处 理 , 用 1/4 节点 代替 
裂 尖 节 点 ,以 实现 该 点 处 的 应 力 奇异 。 


习 1 围绕 裂纹 尖端 的 积分 路 径 
Fig.1 Integral path around crack tip 


C Bl A 8B5 
图 2 1/4 奇异 单元 


Fig.2 1/4 singular element 


2 有 限 元 模型 


采用 随机 顺序 吸附 (RSA ) 算 法 ,利用 脚本 语言 
在 ABAQUS 中 生成 SFRC 模型 。 建 模 采 用 的 假设 : 
中 短 纤维 用 弹性 粱 单元 模拟 ; @ 短 纤维 之 间 不 存在 
交叉 ;@) 纤 维 与 基体 界面 完美 , 短 纤维 均 能 人 均匀 的 
基体 中 ;9 短 纤 维 在 基体 中 的 长 度 统一 。 几 何 模型 
了 =50 mm, 玉 =25 mm ,初始 裂纹 长 度 co =6.25 mm， 
方向 平行 于 * 轴 。 含 有 随机 短 纤维 分 布 (V = 
1.0% ,l=1.5 mm,d =0.012 mm) 的 有 限 元 模型 如 
图 3 所 示 。 

本 研究 采用 三 阶 的 Ogden 本 构 模 型 来 描述 基体 
材料 ,应 变 能 密度 函数 由 3 个 主 拉 伸 比 和 A 表示, 即 


3 
1. 
Ws 人 二 人 
i=1 Qi 


其 中 ,橡胶 基体 材料 参数 :ww = -8.116 MPa,a = 
1.975 , 必 =5.725 MPa,a, =2. 388 ,us =5.467 MPa， 
oa = -4.495.7 ; 短 纤维 材料 参数 :已 = 75 000 MPa， 
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=0.3。 采 用 租 入 单元 (embedded element ) 模拟 短 
纤维 和 橡胶 基体 之 间 的 接触 ,该 方法 使 短 纤 维和 橡 
胶 基体 的 网 格 彼此 独立 ,网 格 形 状 规则 ,在 一 定 程度 
上 延缓 了 数值 计算 不 收敛 的 发 生 。 在 试 样 一 端 施加 
x 和 y 方 向 约束 , 另 一 端 施加 沿 y 向 的 位 移 载荷 用 以 
模拟 单 轴 拉 伸 情 况 。 假 设 裂纹 尖端 前 缘 为 半径 r = 
0.28 mm 的 圆 。 


4 » 
PNT [太公 
~ MA 


图 3 随机 分 布 有 限 元 模型 


Fig.3 Random finite element model 
3 ”结果 与 讨论 


3.1 网 格 尺寸 对 积分 的 影响 


网 格 精 度 对 有 限 元 的 近 代 会 产生 一 定 的 影响 ， 
这 不 可 避免 地 会 造成 误差 ,特别 是 裂纹 尖端 处 ,因此 
有 必要 对 单元 网 格 尺寸 进行 研究 ,以 使 其 对 计算 结 
果 的 影响 最 小 。 经 计算 得 出 , 当 裂 纹 前 缘 区 域 以 外 
的 网 格 尺寸 小 于 0.15 mm 时 ,J 积分 值 趋 于 稳定 ,本 
人 研究 取 0. 10 mm。 通 过 在 r=0.28 mm 裂纹 前 缘 区 
域 设置 不 同 的 路 径 数 目 10、15 、20 .25 、30、35 40、45 
和 50 来 分 析 裂 纹 前 缘 网 格 尺 寸 对 /积分 值 的 影响 ， 
结果 如 图 4 所 示 。 在 裂纹 尖端 前 缘 区 域 以 外 的 网 格 
尺寸 为 0.1 mm 的 情况 下 ,裂纹 尖端 前 缘 区 域 网 格 尺 
才 的 改变 对 7/ 积分 有 一 定 的 影响 , 随 路 径 数目 的 增 
加 , 即 裂 尖 前 缘 区 域 网 格 尺 寸 减 小 ,J 积分 值 增 大 ， 
当 路 径 数目 大 于 25 时 ,J 积分 值 趋 于 稳定 。 后 续 研 
究 中 均 和 采用 统一 的 网 格 划分 方式 (裂纹 尖端 前 缘 路 径 
数目 为 35 ,其 余 网 格 尺寸 为 0.10 mm) 以 确保 计算 结 
果 的 准确 性 。 如 Vj=1.0% ,lj=1.5 mm,d=0.012 mm 
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模型 的 网 格 总 数 为 139 421 , 其 中 橡胶 基体 网 格 总 数 
为 130 871( 裂 尖 前 缘 网 格 数目 4 900) ,纤维 网 格 数 
8 550( 单 根 纤 维 网 格 数 为 10) ,橡胶 基体 采用 平面 应 
力 连 续 单 元 (CPS4R ) ,纤维 采用 线性 混合 梁 单 元 
(BH21H) 。 


3.2 积分 路 径 对 / 积分 的 影响 


通过 裂纹 尖端 前 缘 区 域 网 格 的 细 化 处 理 构 造 的 
35 种 积分 路 径 可 用 于 断裂 参量 积分 的 计算 ,计算 
可 得 积分 与 积分 路 径 的 关系 曲线 如 图 5 所 示 。 随 
着 积分 路 径 位 置 与 裂 兴 距离 的 增 大 ,J 积分 值 趋 于 
一 个 近似 不 变 值 。 此 绪论 与 文献 [13 ] 的 结论 一 致 。 
起 初 数据 的 不 稳定 是 由 于 裂 尖 处 应 力 奇 异性 的 存在 
导致 。 为 减 小 误差 ,应 当 丢 弃 与 此 不 变 值 相差 悬殊 
的 围 线 积 分 值 。 因 此 ,本 研究 后 续 数 据 处 理 均 采用 
积分 路 径 11 ~35(x =0.08 ~0.28 mm) 的 J 积分 平 
均值 作为 最 终 计 算 结 果 , 用 于 后 续 分 析 不 同 细 观 变 
量 对 短 纤 维 填充 橡胶 复合 材料 抗 撕 裂 性 的 影响 。 


6 
. 
2 ， 。 
A ge ey 
RD Se 
后 一 a a 
号 2 2 2 里 
23 
本 
站 22 a /二 1.5 mm, V=1.0% 
9 [=1.5 mm, V=0.5% 
2 4 [=2.0 mm, V=1.0% 
10 20 30 40 50 
裂 尖 网 格 尺寸 
图 4 ”网 格 尺 寸 对 积分 的 影响 
Fig.4 Mesh size vs. J integral 
2.5 古国 国 
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网 YY99999999999999999Y1 
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积分 位 置 /mm 
5 积分 路 径 对 /积分 的 影响 
Fig.$5 Path vs. J integral 


3.3 J 积分 与 拉 伸 载荷 的 变化 关系 
短 纤维 随机 分 布 模型 (V,=0.5% ,lj=1.0 mm， 


Es 
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d =0.012 mm) 在 拉 伸 应 变 s 为 10% 、15% 、20% 、 
30% 40% 和 50% 的 条 件 下 进行 计算 ,得 到 初始 边 
缘 直 裂纹 沿 不 同 扩 展 方向 的 积分。 图 6 为 扩展 方 
向 0 示意 图 。 不 同 拉 伸 载荷 作用 下 的 J-9 曲线 如 图 
7 所 示 ,曲线 为 开口 向 下 的 抛物 线 ,最 大 J 积分 值 位 
于 09=0? 位 置 处 。 随 着 拉 伸 应 变 s 的 增加 ,抛物线 趋 
势 越 明显 ,J 积分 值 逐 渐 增 大 。 最 大 J 积分 值 位 于 
0=0" 位 置 处 ,也 就 是 初始 边缘 直 裂 纹 在 起 初 的 扩展 
方向 可 能 沿 着 初始 裂纹 方向 。 


ty 
1 
1 


一 一 一 


图 6 裂纹 扩展 方向 0 示意 图 
Fig.6 Crack propagation angle 0 


p=0.5%, [=1.0 mm, 4=0.012 mm 
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图 7 不 同 拉 伸 载荷 下 的 J-6 曲线 
Fig.7 Tensile load vs. J-0 curve 


3.4 J 积分 与 短 纤维 体积 含量 \ 长 度 、 直 径 的 关系 


为 研究 短 纤维 体积 含量 、 长 度 和 直径 对 橡胶 复 
合 材料 抗 括 裂 性 的 影响 规律 。 本 研究 分 别 计算 了 短 
纤维 体积 含量 V =0,0.005 ,0. 010 ,0.015 ,0. 020 , 短 
纤维 长 度 j=1.0 mm,1.5 mm,2.0 mm,3.0 mm, 短 
纤维 直径 d = 0. 012 mm , 0. 016 mm,0. 020 mm， 
0.026 mm ,0.030 mm ,0. 036 mm 时 橡胶 复合 材料 沿 
着 扩展 方向 9=0° 的 积分 。 当 短 纤 维 长 度 /= 
1.5 mm, 短 纤维 直径 d 对 积分 的 影响 如 图 8(a) 所 
示 , 可 以 看 出 随 着 短 纤维 直径 的 增加 ,J 积分 值 减 
小 , 即 橡胶 复合 材料 的 抗 撕 裂 性 得 到 提高 。 图 8(b) 
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显示 了 d = 0.012 mm, 沿 初始 裂纹 方向 扩展 的 / 积 
分 值 随 着 短 纤 维 体积 含量 V 和 短 纤维 长 度 二 的 变 
化 ,可 以 看 出 短 纤 维 的 加 入 使 /积分 值 明 显 降低 , 即 
短 纤 维 可 以 提高 橡胶 复合 材料 的 抗 括 裂 性 。 随 着 短 
纤维 体积 含量 的 增 大 或 短 纤维 长 度 的 增加 ,J 积 4 
值 逐 渐 减 小 ,橡胶 复合 材料 的 抗 撕 裂 性 逐渐 增 大 。 

短 纤维 的 加 入 可 以 提高 橡胶 复合 材料 的 抗 折 
性 是 因为 短 纤 维 分 子 中 的 茶 环 较 大 ,其 刚性 和 模 量 
较 大 ,加 之 直径 较 小 、 比 表面 积 较 大 ,在 橡胶 发 生 折 
裂 时 ,可 阻止 裂纹 扩展 ,纤维 从 橡胶 基体 中 拔 出 时 消 
耗 大 量 的 能 量 ,因此 表现 出 较 高 的 撕 裂 强度 。 在 
SFRC 复合 材料 中 ,裂纹 的 扩展 因 短 纤维 的 阻碍 而 改 
变 扩 展 方向 ,而 随 短 纤维 含量 的 增加 ,裂纹 频繁 改变 
扩展 方向 ,导致 扩展 路 径 变 长 ,曲折 度 增 大 ,整体 试 
样 拉 断 时 所 需要 消耗 的 功 比 裂 纹 扩 展 取 直线 路 径 时 
要 大 ,因此 抗 撕 裂 性 得 以 提高 。 随 着 短 纤维 长 度 和 
直径 的 增 大 , 短 纤维 与 橡胶 基体 相互 作用 的 表面 积 
越 大 ,纤维 在 裂纹 扩展 过 程 中 对 橡胶 基体 束缚 作用 
越 强 ,在 橡胶 发 生 撕 裂 时 ,可 有 效 抑 制 裂纹 的 扩展 ， 
因此 表现 出 较 高 的 抗 撕 裂 性 。 
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图 8 J 积分 随 短 纤维 细 观 参量 的 变化 


Fig.8 Microparameters of short fibre vs. J integral 
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3.5 J 积分 与 短 纤维 取向 角 的 变化 关系 


在 SFRC 混 炼 过 程 中 ,橡胶 的 理性 流动 和 炼 胶 
机 辊 简 的 剪 切 拉 伸 强度 效应 , 容易 使 得 刚性 短 纤维 
趋 于 定向 排 布 。 已 有 研究 表明 , 短 纤维 取向 后 的 胶 
料 平行 于 纤维 取向 方向 的 机 械 性 能 远 远 高 于 垂直 于 
纤维 取向 方向 的 ,其 复合 材料 具有 各 向 异性 。 张 立 
群 等 5 通过 实验 法 研究 了 SFRC 的 撕 裂 强度 与 短 纤 
维 用 量 关系 ,证 明了 短 纤维 取向 的 各 向 异性 导致 . 
向 和 7 向 撕 裂 机 理 的 不 同 。 

本 研究 采用 两 种 不 同 的 取向 角 分 布 来 研究 短 纤 
维 取 向 角 a 对 橡胶 复合 材料 抗 据 裂 性 的 作用 规律 。 
取向 角 a 在 某 一 角度 范围 内 随机 分 布 ;@ 取 向 角 
a 服从 高 斯 分 布 (w,o?) , 短 纤维 取向 的 高 斯 分 布 函 
数 可 表示 为 


p(0) = 


1 四 县 
a 
式 中 :y 为 短 纤维 取向 角 的 均值 ;o 为 短 纤维 取向 角 
分 布 的 方差 ,本 研究 取 o =50。 通 过 调整 w 值 实现 
取向 角 的 变化 ,使 短 纤维 趋 于 定向 分 布 。V = 
0. 5% ,l=1.0 mm,d =0.012 mm 两 种 纤维 取向 角 分 
布 的 模型 如 图 9 所 示 , 取 向 角 a 指 短 纤 维 与 y 轴 负 
向 的 夹 角 。 


NA A 人 心 攻 We Hl 
NR AN IN YA | NA 
NO a Pp bd 

CR A 
小 DA 
WA AS 
1 和 i 1、 
AN \ 
VA # #1 | 
RA 
NN WR UN ) " MA A 
N NW A \ / 为 
N Ra 由、 从 WW 
\ RN a1 中 1 
| 了 WE 4 
\ yl I'\ 1 
SN 
和 全 1 
人 
wera 
AN Wy A 
人 
eA i 
\ 4 I YY 
PE NN 
\ 1/ 人 
人 AN 全 
1 A 和 二 1 Vi A 
1 fyy 
WR nl yA { 
4 得， 有 \\ 
\ ha in A A ye 
py 出 是 ， 人 ey 
了 
1 1 1 yy 


(a) 高 斯 分 布 (30°) 


(b) 随机 分 布 (-30°~30°) 
图 9 短 纤 维 定向 分 布 有 限 元 模型 
Fig.9 Aligned fiber model 

短 纤维 取向 角 a 服从 高 斯 分 布 的 /积分 值 随 取 
向 角 的 变化 如 图 10 所 示 。 当 取向 角 均 值 从 -90° 变 
化 至 90" 时 ,7 积分 值 呈 先 减 后 增 的 趋势 ,最 小 值 位 
于 0° 位 置 。 当 j=0°、o =50, 由 高 斯 分 布 的 性 质 
知 ,P1IX-Al<o =Pllal<7 =0.68, 即 取向 角 
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a 有 68% 的 机 会 落 在 ( -7°,7°) 范 围 内 。 而 短 纤维 
取向 角 越 接近 横向 纤维 (a =0°) 的 取向 ,复合 材料 
抗 撕 裂 性 越 好 。 这 是 由 于 当 和 裂纹 扩展 过 到 纵向 纤维 
(a=90°) 时 ,扩展 方向 沿 短 纤维 的 取向 进行 , 绕 行 
距离 短 ,消耗 能 量 较 小 , 抗 撕 裂 性 较 弱 。 而 横向 纤维 
会 阻碍 裂纹 的 扩展 ,使 扩展 路 径 曲折 度 增加 , 绕 行 距 
离 变 长 , 消耗 更 多 的 能 量 , 从 而 提高 材料 抗 撕 裂 
性 。 

短 纤维 取向 角 在 - 10° ~ 10°%、- 30° ~ 30°、 
-50° ~5$0"、-70" ~70" 和 -90" ~90°? 这 5 个 范围 内 
随机 分 布 的 7 积分 值 如 图 11 所 示 , 随 取向 角 角 度 范 
围 的 增 大 ,J 积分 值 逐 渐 增 大 , 即 材料 发 生 撕 裂 的 可 
能 性 逐渐 增 大 。 这 是 由 于 取向 角度 范围 越 小 , 短 纤 
维 越 趋 于 定向 分 布 , 短 纤维 取向 角 与 横向 纤维 (a = 
0°) 夹 角 越 小 ,横向 纤维 数目 越 多 ,而 横向 纤维 对 裂 
纹 扩 展 的 阻碍 作用 更 为 显著 。 随 着 短 纤维 取向 角度 
范围 的 增 大 ,横向 纤维 数量 减少 , 当 取 向 角度 范围 为 
-90° ~90°( 即 完全 随机 取向 ) ,此 时 横向 纤维 数目 
最 少 , 抗 撕 裂 较 弱 。 


3.0 p=0.5%, 7=1.0 mm, d=0.012 mm 
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Fig. 10 Fiber orientation angle vs. J integral 
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Fiber orientation angle vs. J integral( random oriented) 
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3.6 实验 验证 


委托 某 轮 胎 胶 料 生产 三 家 制备 试 样 。 采 用 美国 
杜邦 公司 Kevlar 芳 纶 短 切 纤维 制备 SFRC 试 件 , 基 
本 配方 见 表 1。 对 短 纤维 外 观 形 貌 、 试 件 横 断面 、 纵 
断面 与 厚度 方面 垂直 的 截面 进行 电镜 观测 。 对 电 
镜 图 像 观测 并 分 析 可 以 得 出 以 下 结论 :QD 图 12(a) 
为 Kevlar 纤维 的 外 观 形 貌 ,直径 大 约 为 0.012 mm; 
@ 图 12(b) 和 图 12(c) 中 有 不 少 纤维 从 面 内 冒 出 ， 
而 图 12(d) 呈现 的 平面 较为 光滑 ,基本 无 纤维 从 面 
内 冒 出 ,这 说 明 纤 维基 本 不 治 试 件 的 厚度 方向 取向 ， 
这 是 由 于 在 试 片 成 型 过 程 中 , 短 纤维 刚性 较 大 ,容易 
沿 胶 料 流动 面 内 取向 。 实 验 中 试 样 选用 沿 胶 料 流动 
方向 切割 的 哑铃 型 试 件 (简称 X 试 样 ) ,在 试 件 中 部 
预制 边缘 裂纹 ,采用 由 团队 研发 肯特 机 械 电 子 有 限 
公司 生产 的 电子 拉 伸 机 进行 单 轴 拉 伸 。 


500 um 


世 1 A 


300 hm 


(0) 试 件 纵 断 (qd) 与 厚度 方向 垂直 的 截面 
到 12 短 纤维 外 观 形 貌 及 试 件 内 部 结构 电镜 观测 图 
Fig.12 The SEM of Kevla fiber and the internal 


structure of specimen 
表 1 SFRC 配方 
Tab.1 Ingredient of SFRC 


配方 成 份 (每 百 份 橡胶 含量 ) X 试 样 
天 然 橡 胶 80 
丁 葵 橡胶 -1502 20 
芳 纶 短 切 纤维 所 
中 超 耐 磨 炭 黑 N -220 50 
其 他 配合 剂 18.4 
纤维 长 度 /mm 4 
纤维 体积 分 数 0.0317 


在 文献 [17-19 ] 中 已 经 介绍 了 准 静 态 裂 纹 萌 生 
和 扩展 的 阻力 的 详细 测试 过 程 。 简 而 言 之 ,将 具有 
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不 同 预 切 割 长 度 的 试 样 以 5 mm/min 的 恒定 十 字 头 
速度 拉 伸 。 根 据 下 式 计 算出 在 一 定位 移 下 具有 一 定 
裂纹 长 度 的 试 样 的 值 。 


于， 


其 中 :U 是 在 一 定位 移 A 下 的 应 变 能 ;4 是 断裂 表面 
积 ;8 和 a 分 别 是 厚度 和 预 切 割 长 度 。 

具有 不 同 初始 裂纹 长 度 单 轴 拉 伸 至 断裂 下 相应 
的 载荷 -位 移 曲 线 如 图 13 所 示 。 对 于 任何 给 定 的 裂 
纹 长 度 ,在 断裂 起 始 前 ,通过 积分 相应 的 载荷 -位 移 
曲线 ,在 不 同 的 位 移 值 下 评估 能 量 U( 图 14)。 图 15 
为 在 不 同位 移 水 平 w 下 ,断裂 试验 中 的 输入 能 量 与 
裂纹 长 度 的 关系 ,通过 对 数据 进行 线性 插值 ,可 以 得 
到 一 系列 斜率 为 dU/dal, 的 直线 ,代入 式 (5 ) 即 可 得 
到 J 值 ,实验 结果 列 于 表 2。 
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图 13 单 轴 拉 伸 至 断裂 试验 的 载荷 -位 移 曲 线 


Fig. 13 Load-displacement curve of uniaxial 


tensile to fracture test 


0 1 2 3 4 5 6 7 8 
位 移 /mm 
图 14 不 同位 移 水 平 u 下 的 能 量 评估 


Fig.14 Energy assessment at different displacement levels u, 


由 作者 先前 工作 :加 可知 X 试 样 中 短 纤维 取向 


9 10 1 


角 趋 近 于 高 斯 分 布 。 采 用 j=0° .go? =50 建立 = 
0.0317,4j=4.0 mm,d=0.012 mm 的 高 斯 分 布 模型 
进行 /积分 值 计 算 ,数值 结果 如 表 2 所 示 。 可 见 在 
小 应 变 下 ,实验 值 和 仿真 值 较 接近 ,验证 了 模型 的 可 
靠 性 。 造 成 误差 的 原因 是 数值 模拟 中 未 考虑 界面 缺 
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陷 ,以 及 随 着 纤维 长 度 的 增加 ,由 于 纤维 柔韧 性 更 容 
易 发 生 弯曲 ,在 这 种 情况 下 基于 长 直 纤 维 假设 得 到 
的 数值 结果 会 偏离 实验 结果 。 


140 


100 u/mm 


初始 型 纹 长 度 /nm 
图 15 断裂 过 程 输入 能 量 与 初始 裂纹 长 度 的 关系 
Fig.15 The input energy vs. initial crack length 
表 2 ,/ 积分 值 实验 值 和 仿真 值 对 比 


Tab.2 Comparison of experimental value and 


simulation value of J integral 


位 移 值 /mm 1 2 3 4 
实验 值 /(N . mm-!) 0.45 1.49 2.92 4.88 
仿真 值 /CN . mm-!) 0.43 1.23 2.82 5.56 

相对 误差 0.04 0.17 0.03 0.14 

人 
4 结 论 


本 研究 采用 ABAQUS 软件 计算 了 SFRC 随机 分 
布 模型 沿 不 同 扩展 方向 的 /积分 ,研究 了 网 格 数目 
和 积分 路 径 对 /了 积分 计算 结果 的 影响 ,分 析 了 复合 
材料 抗 据 裂 性 随 短 纤维 含量 长度 .直径 变化 关系 ， 
讨论 了 取向 角 对 SFRC 抗 撕 裂 性 的 影响 以 及 短 纤 维 
在 某 一 角度 范围 内 随机 取向 的 /积分 值 随 角 度 范 转 
的 变化 情况 。 数 值 模拟 得 到 J 积分 值 与 试验 结果 比 
较 接近 ,验证 了 数值 模型 的 合理 性 和 可 靠 性 。 得 出 
如 下 结论 。 

1) 裂纹 前 缘 区 域 网 格 尺 寸 和 积分 路 径 对 了 积分 
的 数值 计算 结果 有 一 定 的 影响 。 随 着 网 格 尺 二 的 减 
小 和 路 径 位 置 与 裂 尖 距离 的 增 大 ,7 积分 值 趋 于 稳 
定 。 为 了 确保 结果 的 精度 和 可 比 性 ,网 格 划 分 方式 
应 一 致 . 网 格 尺 二 应 足够 小 .积分 路 径 应 远离 裂纹 尖 
端 以 避免 其 对 积分 值 计算 的 影响 。 

2) 同 一 边缘 直 和 裂纹 (ao =6. 25 mm ) 的 / 积分 随 
拉 伸 载荷 的 增 大 而 增 大 ,在 相同 拉 伸 应 变 下 , 沿 初始 
裂纹 方向 的 了 积分 值 随 短 纤维 含量 .长 度 和 直径 的 
增 大 而 减 小 , 即 橡胶 复合 材料 抗 撕 裂 性 能 随 短 纤维 
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含量 .长度 和 直径 的 增 大 而 提高 。 

3) 短 纤维 在 橡胶 基体 中 的 取向 分 布 会 影响 橡胶 
复合 材料 的 抗 据 裂 性 。 取 向 角 a 服从 高 斯 分 布 (x， 
0 ) 人 情况 下 ,J 积分 值 随 取 向 角 期 望 值 的 增 大 而 增 
大 。 当 取向 角 均 值 为 0° ,J 积分 值 最 小 。 这 是 由 于 
取向 角 为 0" 的 纤维 ( 横向 纤维 ) 会 阻碍 裂纹 的 扩展 ， 
使 得 绕 行 路 径 曲折 , 而 绕 行 距离 越 长 , 消耗 能 量 越 
多 , 抗 撕 裂 性 越 好 。 

4) 当 短 纤维 在 一 定 角度 范围 内 随机 取向 ,取向 
越 大 , 短 纤 维 取向 分 布 越 随机 。, 积分 值 随 着 
可 范围 的 的 增 大 而 增 大 , 即 短 纤维 在 小 角度 范围 
机 取向 ( 趋 于 定向 取向 ) 的 抗 扬 裂 性 优 于 大 角 
站 内 随机 取向 ( 趋 于 随机 取向 )。 
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